Relevance of coordinate and particle-number scaling in density functional theory 
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We discuss a /3-dependent family of electronic density scalings of the form n\(r) = A 3/3+1 n(A^r) 
in the context of density functional theory. In particular, we consider the following special cases: the 
Thomas- Fermi scaling (/3 = 1/3 and A ^> 1), which is crucial for the semiclassical theory of neutral 
atoms; the uniform-electron-gas scaling (/3 = —1/3 and A > 1), that is important in the semiclassical 
theory of metallic clusters; the homogeneous density scaling (/} = 0) which can be related to the 
self-interaction problem in density functional theory when A < 1; the fractional scaling (/? = 1 and 
A < 1), that is important for atom and molecule fragmentation; and the strong-correlation scaling 
(P = — 1 and A > 1) that is important to describe the strong correlation limit. 

The results of our work provide evidence for the importance of this family of scalings in semi- 
classical and quantum theory of electronic systems, and indicate that these scaling properties must 
be considered as important constraints in the construction of new approximate density functionals. 
We also show, using the uniform-electron-gas scaling, that the curvature energy of metallic clusters 
is related to the second-order gradient expansion of kinetic and exchange-correlation energies. 

PACS numbers: 71.10.Ca,31.15.E-,31.10.+z 



I. INTRODUCTION 

In Density Functional Theory (DFT) Q-ij], the density 
scaling is a key concept and was used along the years to 
derive many exact constraints P4l2|. as well as useful 
virial relations [J, for the kinetic, exchange and 

correlation energy functionals. Furthermore, the density 
scaling has a fundamental role in DFT due to its intimate 
relation with the adiabatic connection formalism 0, EH, 

HJ. 

One important family of scaling transformations for 
the particle density n is defined by the general linear 
transformation [5j] 



n M , a (r) = det(M) n(Mr - 



(1) 



where M is a real invertiblc 3x3 matrix and a 6 R 3 . Scal- 
ing transformations of this kind correspond to changing 
the external potential associated with the density whilst 
preserving the normalization of the density, i.e. 



J n M , a (r)dr = J n(r)dr = N 



(2) 



The most important scalings connected with Eq. ([1} arc 
the ones defined by 



Mtj = A , Sij , Xj > 



and a = 



(3) 



i.e., with Mij a diagonal matrix with positive elements. 
In particular, there are three cases of high physical inter- 
est: 

(i) the uniform scaling defined by Ai = A2 = A3 = 
A > 0, hence n\(r) = X 3 n(Xr). Under uniform scaling 
the Kohn-Sham exchange and the non-interacting kinetic 
energies transform as E x [n\] = XE x [n) and T s [n\] = 
A 2 T s [n], respectively However, the interacting kinetic 



energy T[n] and the non-local Hartree-Fock exchange en- 
ergy E X F [n] do not have these properties [1| . All popu- 
lar exchange and kinetic energy functionals are nowadays 
designed to satisfy the uniform scaling relation 0, Q . 

(ii) The two-dimensional nonuniform scaling of the 
density defined by Ai = 1 and A2 = A3 = A > 0, hence 
n z ^{x, y, z) = X 2 n{x, Xy, Xz), under which the system ap- 
proaches the ID limit when A — > 00. Although there are 
important conditions for both Kohn-Sham exchange and 
correlation energies under this scaling [5|, they are not 
satisfied by any popular XC functional (to our knowl- 
edge). 

(iii) The one- dimensional nonuniform scaling of the 
density defined by Ai = A2 = 1 and A3 = A > 0, 
hence n\(x,y,z) = Xn(x,y, Xz), under which the sys- 
tem approaches the 2D limit when A — > 00. This 
density scaling, that is related to the dimensional 
crossover of the XC energy (from 3D to 2D) (H^U, 
has been recently incorporated in a semilocal XC func- 
tional, named q2D-generalized-gradient-approximation 
(q2D-GGA) [25[, constructed for mild and strong quasi- 
2D regimes. The remarkable performance of the q2D- 
GGA for surface energies and lattice constants of tran- 
sition metals, showed the power of the one-dimensional 
nonuniform scaling. 

In this paper we consider a different type of scaling 
for the density and focus the attention on the family of 
scaling relations of the form 



n x (r) = A 3/?+1 ? i(A^r), A > 0, 



(4) 



where j3 is a parameter. These scaling transformations 
differ from those defined by Eq. (|T|) in the fact that 
they do not only change the external potential associ- 
ated with the density n, but also provide a change in the 
particle number (N —5- XN). Well known members of 
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the scalin g fa mily defined by Eq. ((4]) are the Thomas- 
Fermi [2II l27j and the homogeneous density scaling 
which are relevant for the semiclassical theory of the 
many-electron, non-relativistic, neutral atom [26[ , and 
for the self- interaction error [9(, respectively. (Unless 
otherwise stated, atomic units are used throughout, i.e., 
e 2 = h = m e = 1.) 

This article is organized as follow: in Section [TT| we 
shortly review the fractional-particle density functional 
theory, that is the correct framework for the scalings of 
Eq. (Q}; in Section IIIII we present scalings properties 
of useful density functionals; Section IIVI is devoted to 
the physical properties of selected density scalings (for 
,3=1/3, -1/3, 0, 1, and -1); and in Section |V] we briefly 
analyze the performance of popular density functionals 
for the above scalings. Finally, in Section l*VTl we summa- 
rize our conclusions. 



II. THEORY AND DEFINITIONS 

Under the scaling transformations defined by Eq. Q 
the normalization of the particle density is modified, 
thus the number of particles in the system is varied as 
N — > XN, with A-/V being in general a non-integer num- 
ber. The conventional picture of DFT, based on the 
Hohemberg-Kohn theorems [28[ and/or the Levy con- 
strained search 29] is thus inappropriate in this case. In 
fact, a proper definition of densities with a non- integer 
number of particles is required in this case [30l - [3^ | . In this 
work we achieve this through the introduction of ensem- 
bles densities within a zero-temperature grand canonical 
ensemble theory [30l |. 

The central quantity to consider is therefore the 
density-matrix operator [3[ 



Mi 



(5) 



where the sums are extended over all possible particle 
numbers M and over all the states of the M-particle 
Hamiltonian, \^>Mi) is the i-th eigenstate of the M- 
particle Hamiltonian and pui is the probability weight 
to find the system in the eigenstate \^>Mi) (hence, we 
have < pm% _ 1 and Y^MiPMi — !)• The expecta- 
tion value of any operator A is consequently obtained as 
Tr TA . In particular, for the particle density we have 



n(r) = Tr fn(r) = ^p fe n fe (r) , 



(6) 



where we defined the super-index k — Mi, n k is the pure- 
state density of the M-particle z-th state, and we used 



i(r) = f/S + (r)^(r) , $(r) = ^<^-(r)c 



(7) 



with 4>i and a, the single-particle orbital and annihilation 
operator of state i, respectively. The total number of 



particles is, according to Eq. (J6)), 

N = Y,PkN k , N k =Jn k (r)dr. (8) 

k 

Similarly, for the one-particle density-matrix operator we 
have 7i(ri,r 2 ) = ^ (ri)^(r 2 ) and 



7i(ri,r 2 ) = Tr Fyifr^ra 



X^fc7ifc(ri,r 2 ) ■ (9) 

k 



The ensemble Hohemberg-Kohn universal functional 
and the non-interacting kinetic energy are defined as [|[ 



F[n] 



min Tr 



T s [n] = min Tr 
f-m 



r 

TT 



(f + U ee ) 



(10) 
(11) 



where V ee is the electron-electron repulsion operator and 
the kinetic energy operator is defined by 

f = y [V? a 7i(ri,r 2 )](S(ri-r 2 )dndr 2 . (12) 

The electron-electron repulsion operator can be further 
decomposed into Coulomb, exchange, and correlation 
contributions to yield 



J 



1 f n(ri)n(r 2 ) 



E x [n] = Tr 
E c [n] = Tr 



ri -r 2 



r minE c 



dri<ir 2 , 



(13) 
(14) 
(15) 



where T m i n is the density-matrix operator minimizing 
either f + V ee or f, according to Eqs. (dHJ) and (fTTj) , and 
with the exchange energy operator defined as [3|] 



E , 



|-7i(ri, r 2 )|- 
ki -r 2 | 



dridr 2 



(16) 



while no explicit expression is known for E c . Note that 
all the traces can be easily evaluated by use of Eqs. 
^ and ©, together with the resolution of identity 

I = 5lMi\*Mi)(*m\. 

Within the theoretical framework sketched above the 
scaling relations of Eq. (Q| can be interpreted as an 
uniform scaling of the pure-state densities (n^(r) — > 
A 3/3 rife(A' 3 r) Vfc) accompanied by a remodulation of the 
statistical weights such that the particle number is 
changed to A7V. Note that this latter is in general 
a complicated transformation because it must accom- 
plish the required particle-number variation preserving 
the correct density and without violating the normaliza- 
tion conditions for the statistical weights (0 < pk < 1 
and J2kP k = !)■ l n particular, the simple transforma- 
tion pk — ► Xpk in general is not a suitable transforma- 
tion, as it brings a violation of the normalization condi- 
tions. Nevertheless, one special case is when we con- 
sider a system with one electron or less and A < 1. 
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In this case in fact the density of fractional charge q 
can be written n q {r) = qn±(r) + (1 — q)n (r), where 
n\ and hq are the densities for one particle and zero 
particles, respectively (the latter is of course identically 
zero everywhere) . The scaling transformation yields then 
n qX = \q\ 3 P ni (\Pr) + (1 - Xq)X^n (X l3 r). Thus, in this 
special case we have indeed p — > Xp (all information con- 
cerning the no term can be neglected). This important 
result will be employed in next section to derive exact 
scaling relations for the non-interacting kinetic energy 
and the Kohn-Sham exchange for systems with fractional 
occupation. 

More insight into the properties of the transformation 
governing the statistical weights can be achieved by con- 
sidering (in analogy with the usual adiabatic connection 
procedure) the M-particle Hamiltonian with a local po- 
tential which gives the right pure-state densities to re- 
cover Eq. ©. In this way the trasformation of the sta- 
tistical weights, connected to scaling transformations of 
the type denned in Eq. Q, can be defined explicitly 
in two separate cases. When, upon scaling, the vari- 
ation in the particle number is smaller than one (i.e., 
1 < A< (N + l)/N or similarly (N - l)/N < X < 1), the 
statistical weights change with A as 

p N : 1 1 - N(X - 1) ; p N+1 :0^N{X- 1), (17) 

and consequently, the exact total energy functional E[n] 
varies linearly with A j30j 

E[n x ] = (1 - 7V(A - 1)) E[n N ]+N{\-l)E[n N +i] , (18) 

where njv and tijv+i denote N- and (JV + l)-densities. 
When the variation in the particle number is larger than 
one (i.e., the number of particles changes from N to N + 
L + uj, with L an integer and \u\ < 1), the statistical 
weights change as 

Pn ■ 1 ; p N+L :0^1-uj ; Pn+l+i ■ -> w . 

(19) 

Therefore, the transformation can be seen as a chain 
of transformations like the one in Eq. concern- 
ing successively p N , pn+i, Pn+l, Pn+l+i- As a 
consequence, the exact total energy functional will be 
described as a succession of straight lines. Thus, the 
derivative discontinuity [30] plays an important role for 
the scaling transformations, especially when TV and L are 
finite. However, in the limit L — > oo the role of derivative 
discontinuity is diminished and one can always consider 
u) = 0. In this case the behavior of the system is closely 
related to the semiclassical physics, as shown in the next 
sections. Because of the complexity of this density scal- 
ing, in this paper we mainly consider only two extreme 
cases of interest: when N = 1 and A < 1, and the semi- 
classical limit A -+ oo. 

III. GENERAL SCALING PROPERTIES 

Using the definition of Eq. (Q} we can immediately 
find the scaling properties of several quantities which de- 



pend explicitly on the particle density. This is the case, 
for example, of the the Coulomb energy and any exter- 
nal potential (e.g., E nuc = J nv nuc dr with v nuc oc 1/r). 
Hence, 

J[n x ] = A^+ 2 J[n], (20) 
E nuc [n x ] = X (3+1 E nuc [n]- (21) 

Another interesting quantity is the local Seitz param- 
eter r s = [3/(47rn)] 1 / 3 which scales as 

r sX (r) = X-P-h s (XPr) . (22) 

Therefore, for j3 > —1/3 we have that A — > oo implies 
r s \ — > 0, while A — > implies r s \ — » oo. Thus, the 
conditions A — 5- oo and A <C 1 correspond to scalings 
to the high- and low-density limits, respectively. The 
opposite is true for (3 < —1/3. For the special case (3 — 
1/3 instead the local Seitz parameter is independent on 
A and the density regime cannot be modified by a scaling 
transformation. 

On the other hand, the usual density parameters s — 
\Vn\/2k F n, q = V 2 n/{4(37r 2 ) 2 / 3 n 5 / 3 }, t = \Vn\/2k s n, 
and v = |Vn|/2fc„n, with k F = (37r 2 n) 1 / 3 being the lo- 
cal Fermi wave- vector [3^, [34[, k s = [Akphr) 1 /' 2 being 
the Thomas- Fermi screening wave- vector [33|, l34j . and 
k v = 2(3/(47r 4 )) 1 / 18 n 1 / 9 being the wave vector suitable 
for bonding and valence regions (35j, scale according to 

8 X (t) = X~h (X r) , qx(r) = X-iq(X r) , (23) 

t x (r) = A^-h(A^r) , v x {r) = A^"*u (A"r) (24) 

Thus, the reduced gradient and Laplacian for exchange 
and kinetic energies (s and q) are independent on j3 so 
that the slowly-varying density limit (s, q — > 0) is reached 
whenever A — > oo, while for A <C 1 a rapidly- varying 
density regime is always set up. On the contrary the 
density parameters t and v, which are relevant for the 
correlation, have a dependence on ft. Therefore, they 
can describe different density regimes depending on the 
actual value of the parameter j3. 

Because for any value of (5 the reduced gradient for ex- 
change and kinetic energy s and the reduced Laplacian 
q become small in the limit A — > oo, it is also interest- 
ing to investigate the scaling behavior of local density 
approximations (LDA) and gradient expansions for the 
non-interacting kinetic energy and the exchange energy. 
These expressions will in fact become almost exact in the 
limit A — > oo. The required scaling relations are given by 



the formulas 






E L x DA [nx] 


= X^+ 4 / 3 E^ DA [n], 


(25) 


E° E2 [nx] 


= X^ 3 E^ E2 [n], 


(26) 


T s LDA [nx] 


= X 2 ^ 3 T s LDA [n], 


(27) 


T? E2 [nx] 


= (l/9)Tf [n x ] = {l/9)X 2fs+1 T^ 


[n],(28) 


T? Ei [nx] 


= A^ +1 / 3 T s G£4 [n] , 


(29) 



where E^ E2 is the second-order gradient correction 
(GE2) term of the exchange energy [36], Tf E2 is 
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the second-order gradient correction term of the non- 
interacting kinetic energy Isj, is the von Weizsacker 
kinetic energy functional [2], and T^ E4 is the fourth- 
order kinetic energy gradient expansion term (GE4) 
|38| . Using these expressions, as well as Eq. (|20[) . it is 
possible to provide an useful accurate approximation for 
the universal functional of Eq. (|10p in the slowly- varying 
limit (A — > oo): 

F[n x ] w X^ +2 J[n}+X 2fs+5 / 3 T s LDA [n} + 

+ X 2 P +l T GE2 [n] + X 2 P +l/3 T GEi [n] + 

+X^ 3 E^ DA [n} + X^ +2 / 3 E^ E2 [n] . (30) 

In this formula we neglected correlation contributions. 
In fact, the leading term of LDA correlation energy in 
the high-density limit is ^ E LDA [n] cx J drnln(r s ), 
and does not respect any simple scaling; whereas the 
leading term in the low density limit is [39( E LDA [n] cx 
-/ drnrj 1 , and scales as E LDA [n x ] = X f3+i / 3 E LDA [n). 
Moreover, the second-order correction to the correlation 
energy [H| scales always as Ef m [n\] = X p+2 / 3 E^ E2 [n}. 
However, due to the dependence on /3 of the reduced gra- 
dients t and v, the slowly- varying limit is only reached for 
< 1/6. In addition, for any j3 > —1/3 the scaling to the 
slowly-varying limit corresponds also to a scaling to the 
high-density limit (r sX — > when A — > oo). In this limit 
the correlation contributions are negligible with respect 
to the exchange part, and thus Eq. ([ST)]) becomes almost 
exact. For f3 < —1/3 however correlation corrections to 
Eq. ([30]) might be needed. 

To conclude this section we consider briefly the cases 
of Kohn-Sham kinetic and exchange energies. Using the 
convexity arguments derived in Ref. and taking into 
account that the scaling family of Eq. (QJ can be seen 
as an uniform scaling followed by a homogeneous scaling 
(i.e. n\(r) = X(X 3 °n(X°r))), the following inequalities 
hold: 

T a [nx] > X 2f)+1 T s [n], A > 1, (31) 
T s [n\] < X 2fi+1 T s [n], A < 1, (32) 

and 

\E x [n x }\ > X p+1 \E x [n]\, A > 1, (33) 
\E x [nx]\<X p+1 \E x [n]\, A < 1. (34) 

Moreover, using the rigorous bound T s > Tf, and the 
one conjectured by Lieb H3 T s < T EDA + Tf (for a rig- 
orous, and tighter upper bound of T s in terms of T^ DA , 
Tf and TV, see Eq. (23) of Ref. 0), we can easily 
derive the following inequalities 

X ^+l T W [n] < Ts[nx] < A 2^5/3 T LDA [n] + X 2P + l T W [n] ^ 

(35) 

for any A and (3. In the case of exchange energy, the 
Lieb-Oxford bound HifH E x > E xc > 2.27E LDA gives 



In the special case of one particle or less {N < 1) and 
A < 1, using the formalism presented in Section [II] (es- 
pecially the fact that p — > Xp) it can be shown (see Ap- 
pendix^ that 

T s [n x ] = X 2 ? +1 T s [n], (37) 

and 

E x [n x ] = X^ +2 E x [n}. (38) 

These results correctly agree with the scalings of J[n] 
and Tj^[n], because in case AiV < 1, E x [n] = — J[n] and 
T s [n] =_T] v [n] 53 . (Note also that Eq. (JSTJ) holds for 
N < 2 [if .) 

IV. SELECTED SCALING RELATIONS 

In this section we analyze in more detail the scaling 
relations corresponding to special values of the parameter 
13. In this way we can highlight the physical significance 
of the family of scaling transformations defined by Eq. 
((H) and provide evidence for its importance in electronic 
structure theory. 

A. Thomas- Fermi scaling (/3 = 1/3) 

If we require the LDA non-interacting kinetic energy to 
scale as the Coulomb energy, we find f3 — 1/3, which cor- 
responds to the well known Thomas- Fermi scaling [26| . 
With this choice the scaling of the density and the re- 
duced gradients are 

n x (r) = A 2 n(A 1/3 r) , r sA (r) = A- 2 / 3 r s (A x / 3 r](39) 
SA (r) = A- 1 / 3 .s(A 1 / 3 r) , g A (r) = A^^A^r) ,(40) 
t A (r) = i(A 1/3 r) , v\(t) = X 1/9 v(X 1/3 r) . (41) 

Thus, for A — > oo the high-density limit is reached and 
the exchange and kinetic energies, whose behavior is con- 
trolled by the density parameters s and q, are in a slowly- 
varying density regime. Therefore, the universal func- 
tional can be written [26| 

F[n x ] w A 7 / 3 (T S LD >] + J[n)) + 

+x m(rGE2 [n]+E LDA [n]) + 

+X(T^[n]+E^ E2 [n}) , (42) 

and the Thomas-Fermi kinetic energy (T^ DA ) and the 
classical Coulomb energy are the leading terms in the 
total electronic energy (E nuc scales as A 4 / 3 ). This result 
is very important. In fact, the semiclassical theory of the 
many-electron neutral atom [46T - [48j is based on Eq. (|42[) . 
which leads to the semiclassical asymptotic expansion for 
the kinetic energy 



E x [n x ] > 2.27X 0+4 ^ E^ DA [n}. 



(36) 



T s = c N 7/3 + Cl N 2 + c 2 N 5/3 + ... ; N cc X (43) 



5 



and a similar one for exchange [49J , that are very accurate 
(ty pica l error of order 0.5% - 0.2% even for real atoms [H, 
49. 50]). Note that the second term in the kinetic energy 
expansion (ciA 2 ), can not be captured by the Thomas- 
Fermi scaling, being a quantum correction. Recently, it 
has been demonstrated that these as ympto tic expansions 
are also important tools in DFT (26l. [5ll— f53j] . as they have 
been used to construct accurate non-empirical exchange- 
correlation [5l[ and kinetic [52| energy functionals. 

Extensions of Eq. (l43l) have been also proposed for 
general ions and atoms [541457] , which constitute a more 
challenging problem than neutral atoms alone. However, 
even for the first ionization potential of many-electron 
atoms, the extented semiclassical Thomas-Fermi theory 
shows serious drawbacks and limitations (58| . and ac- 
curate results can be obtained only within Kohn-Sham 
DFT or other orbital-dependent schemes. 

Concerning correlation, we can obtain some insight by 
considering that, any reasonable generalized gradient cor- 
rection to the LDA correlation must be designed to can- 
cel the logarithmic divergence of the LDA term under 
uniform scaling to the high density limit [59]. Thus, in 
the high-density limit ef GA cx ln(i 2 ). However, under 
the Thomas-Fermi scaling , while v\ — A 1 / 9 v, so that 
correctly v — > oo for A — > oo (v is a density parameter 
suitable for valence and tail regions that are evanescent 
in a many-electron neutral atom), the density parameter 
t is just independent on A. Hence, the gradient correc- 
tions to the correlation are independent on the scaling 
and the whole correlation energy is dominated by the 
LDA contribution e^ DA oc ln(r s ) for A — > oo. As a con- 
sequence, popular GGA functionals (e.g. PBE [5!|) that 
recover LDA correlation in this limit, can be argued to 
be accurate (exact) [26| . 

Finally, we mention the importance of the Thomas- 
Fermi scaling also for the atomic densities. In fact, under 
a Thomas-Fermi scaling to the high-density limit, the hy- 
drogenic density n = exp(— 2r)/n resembles features of 
the Thomas- Fermi density [60| , becoming slowly- varying 
over a Fermi wavelength (but not over the screening 
length 2ir/k s ). We recall that the Thomas-Fermi den- 
sity, even if does not show shell structure and does not 
decay correctly, is a very good model for the densities of 
heavy atoms [50l . lolTtolij . For an excellent discussion, see 
Ref. m 



B. Uniform-electron-gas scaling (/3 = —1/3) 



reduced gradients behave as 

n\(r) = n(A _1 ^ 3 r) , r sX (r) = r s (\-V 3 r) (44) 
s A (r) = A- 1/3 s(A- 1/3 r) , g A (r) = X^^qiX^^v}^) 
t A (r) = A- 1/3 i(A~ 1/3 r) , v x (r) = \- 1/3 v(\- 1/3 r§46) 

Thus, both the density and the local Seitz parameter are 
independent on A (except for a coordinate scaling), so 
that the scaling does not involve any transformation to- 
wards the high- or low-density limit. On the contrary, 
by construction, all the density parameters vanish in the 
limit A — > oo, so the slowly- varying density limit is fully 
recovered in this case (hence, the name uniform-electron- 
gas scaling). Interestingly, Eqs. (|4"5"j) and (PToT) also show 
that for first-order density parameters (i.e. those depend- 
ing on Vrt) exactly the same dependence on A -1 / 3 is 
found under the uniform-electron-gas scaling. 

For large values of A the universal functional is well 
approximated by 



F[n x ] 



X 5/3 J[n] + 
+\{Tt DA [n] 

+A l/3 (T GB2 [n] . 
+A -1/3 T GB4 W 



E L x BA [n\ 
\+E? E2 \n] 



E. 



LDA 



E, 



n}) + 



(47) 



where the full LDA correlation energy is considered in Eq. 
(|4"T)) . Of course, in the limit A ^ oo the LDA approx- 
imation of F becomes exact (all the density parameters 
vanish in this limit, by construction). In particular, we 
have 



lim E xc [n\] 



\n c 



(48) 



where rtoo(r) = n(0) = const. Note that while this con- 
straint is satisfied by most non-empirical XC semilocal 
functionals (that recover the LDA for a constant density) , 
it can be out of reach for some wavefunction methods, as 
the second-order perturbation theory of M0ller-Plesset 
(MP2) and its modifications 64, 65|, or the random phase 
approximation (RPA) [66J, |67(. However, the sophis- 
ticated orbital-based inhomogeneous Singwi-Tosi-Land- 

Sjolander (ISTLS) method [HlH, as wel1 the XC kernel 
of linear response time-dependent DFT (in the context of 
the adiabatic-connection fluctuation-dissipation theorem 
[HI, [7(J [7l| ) of Ref. 72 1, are accurate for the uniform- 
electron-gas scaling. 

To provide an example of the utility of the uniform- 
electron-gas scaling let us consider neutral jellium clus- 
ters with N electrons and radius R — t^A 1 / 3 , having the 
external potential 



Consider a generic density parameter d oc \Vn\/n a . 
We define the uniform-electron-gas (UEG) scaling, as the 
scaling belonging to the family of Eq. ((4]) that makes d 
small for any value of a in the limit A — > oo. It is easy 
to prove that this scaling is defined by the parameter 
/3 = —1/3. Under such a scaling the density and the 



-Ai, 



' 25*)> r < R 
r>R, 



due to a positive background density 



ra+ (r) = 



_ / 3/47T7- 3 



1° 



r <R 
r > R 



(49) 



(50) 
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FIG. 1: (Color online) Aixr 2 n\ versus the radial distance r for 
the 58e~ Na jellium cluster (A = 1 and A = 92/58), and for the 
92e~ Na jellium cluster (A = l).The areas under the curves 
are the total number of electrons (58 and 92 respectively). 



This external potential has no singularities, so the re- 
duced gradients and Laplacian of the density, are finite 
everywhere inside the bulk. Moreover, the values of 
the density parameters decrease for increasing number of 
electrons N. Indeed, even for intermediate values of N, 
the density is slowly varying over a Fermi wavelength, so 
the extensions of Thomas-Fermi theory become accurate 

[zMEIIl- 

Jellium clusters with different numbers of electrons 
may be thought therefore to be well described by the 
uniform-electron-gas scaling, since under the scaling pro- 
cedure the number of electrons is changed to XN, the 
local Seitz parameter is kept fixed to r s , and the reduced 
gradients are decreased as A -1 / 3 , in full analogy to what 
happens in the jellium clusters. The relation between 
the uniform-electron-gas scaling and the jellium clusters 
can be in fact clearly recognized by a detailed analysis 
of the cluster's electron densities. To this end, in Fig. Q] 
we show the densities of 58e~ and 92e~ Na jellium clus- 
ters, together with the UEG-scaled density of the 58e~ 
cluster where a value A = 92/58 was used. Remarkably, 
the scaled density agrees very well with the 92e~ cluster 
density showing that the variations of the physical prop- 
erties of the different clusters with N can be well captured 
by the uniform-electron-gas scaling. Note that the small 
differences between the scaled and the true density are 
due to quantum oscillations [H[ that are not accounted 
by the scaling procedure, but are of course included in 
the self-consistent Kohn-Sham scheme. These are how- 
ever not very relevant for the analysis that we consider 
henceforth. 

The total energy of a jellium cluster can be written 



E[n] = F[n]+E c 



ct [n} = F[n} + J Vi:l{v)n{v)dv . (51) 



Therefore, using the uniform-electron-gas scaling for A — > 



oo we can write 



E[n] 



X(Tr A [n]+E^ A [n]+E^ A [n)) + 



_ A i/3 (T; 
-0(X 



GE2 r 



E 



GE2\ 



E 



GE2\ 



-1/3^ 



(52) 



where we used the fact the Coulomb and external poten- 
tials cancel each other in the limit of a large number of 
electrons. This expression agrees well with the asymp- 
totic expansion of the total energy of the clusters derived 
from the liquid drop model 73|-75j (recall that N ~ A) 



Airr 3 

E = a—^N + a4irr 2 s N 2/3 + ~f2irr s N 1/3 , 



(53) 



where a is the volume (bulk) energy and a and 7 are the 
surface and curvature energies. A comparison of the first 
terms (those scaling as A) in Eqs. (|52|) and (|53[) shows in 
fact that the uniform-electron-gas scaling correctly yields 
a ~ T S LDA + E^ DA + E^ DA . Moreover, similar with the 
previous Thomas- Fermi scaling (see Eq. (H2])), the sec- 
ond term in Eq. (1531) is a quantum oscillation term, that 
apparently can not be described by the simple UEG scal- 
ing, without a careful analysis of the Friedel oscillations 
near the surface of the cluster. Finally, the third term is 
related to the curvature energy 
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T. 



GE2 



E 



GE2 



E, 



GE2 



(54) 



This is an important result, because the real edge gas 
contains curvature corrections (7o| . that until now had 
not been addressed. 



C. Homogeneous density scaling (/? = 0) 

The homogeneous density scaling [t| [TJ I771 - I82T ] is ob- 
tained from Eq. ((4]) by setting the parameter j3 = 0. 
Under this condition the scaling relations for the density 
and the various density parameters are 



n\(r) = Xn(r) 
s A (r) = X~is (r) 
t x {r) = X-h(r) 



r sX (r) = X-^ 3 r s (r) (55) 
gA(r)=A-fg(r) , (56) 
v x (r)=X- 1 ^v(v) . (57) 



For A — > 00 the high-density slowly-varying limit is ob- 
tained and the universal functional is well approximated 
as 



F[n x ] 



X 2 J[n] + X 5/3 T s L --[n\ 

+ X 4/3 E LDA [n] + XT GE2 [n] _ 
+A 2/3 £ ,GB2 [n]+A l/3 r GB4 [f 



(58) 



In this case the functional of Eq. (|58[) is dominated by the 
classical Coulomb term, and the second leading term is 
the Thomas-Fermi kinetic energy. (Note that the nuclear 
energy grows only linearly with A.) 



More importantly, the homogeneous scaling is a valu- 
able tool to investigate DFT, when the opposite limit, 
i.e. with A < 1, is considered. For this case in fact 
several studies exist on the scaling properties and exact 
constraints of the kinetic and exchange energy function- 
al [H, [ll | a s well as on the static correlation treatment 
in DFT 13, [111 • I n this work we focus instead on the 
role of the homogeneous scaling in the determination of 
the derealization error [85| of the exchange(-correlation) 
functionals. To this end we consider the simple H^~ dis- 
sociation problem, that is associated with the hydrogen 
atom with fractional charge f&ij , and compute 

AE H {q) =E H - E m - E Hl x-o , (59) 

where Eh = E xc + J for the hydrogen atom and < q < 
1 is the partial electronic charge. The quantity AEh 
represents thus the Coulomb and XC energy difference 
between the dissociation of into one hydrogen atom 
plus one proton and that of the dissociation into two 
hydrogen atoms with fractional electron charge q and 1 — 
q. For the exact exchange- (correlation) functional it shall 
be zero at any value of q. However, due to the one- 
electron self- interaction error Q, for approximated XC 
functionals AEh > for any < q < 1, indicating that 
a fractional dissociation is favorable with respect to the 
exact one. 

For a generic exchange functional we can assume (see 
Section© the scaling ££ n3/ [n A ] = \ a E% n y[n], with a < 2. 
Thus, Eq. (l59l) becomes 

AE H (q) = (l-q 2 -(l-q) 2 )J[n H } + 

+ (l-q a -(l-q) a )ET y [n H ] , (60) 

with tih the density of the hydrogen atom and J[njy] = 
0.3125 Ha. Note that in this case correlation plays no 
role, since for q < 1 the correlation energy is zero. Nev- 
ertheless for semilocal functionals, which are not self- 
correlation free, it is also possible to consider an E x ™ y 



TABLE I: Self-interaction error (SIE) as defined in Eq. ([61]) 
for different popular exchange and exchange-correlation func- 
tionals. All values are computed for the hydrogen density and 
expressed in mHartree. 



X-only functionals 

Functional Ref. SIE 

LDAx [86, 87] 44.5 

PBEx [59] 6.6 

APBEx [51] 1.8 

revPBEx [88] 2.0 

PBEsolx [89] 19.8 

PBEintx [90] 16.9 

B88 [91] 2.7 



XC functionals 



Functional 


Ref. 


SIE 


LDA 


[86, 87, 92] 


22.3 


PBE 


[59] 


5.8 


APBE 


[51] 


-3.4 


revPBE 


[88] 


-4.0 


PBEsol 


[89] 


11.9 


zPBEsol 


[89, 93] 


0.0 


PBEint 


[90] 


9.7 


zPBEint 


[90, 93] 


0.1 




FIG. 2: (Color online) Plot of f a (q) = 1 - q a - (1 - q) a for 
several values of the parameter a. 

in place of E x ny . Now, in analogy with Ref. §1, we can 
define the self-interaction error as 

SIE[n] = J[n] +E xc [n] . (61) 

Equation (|60|) then becomes 

AEh (q) = (/a(«) - /<*(<?)) J[n H ] + f a {q)SIE[n H ] , (62) 

with f a (q) = 1 — q a — (1 — q) a . A plot of f a for several 
values of a is provided in Fig. [5] 

Equation (l6"2l shows that AEh is always positive 
whenever the self-interaction error is not zero and the 
scaling behavior of the functional differs from the exact 
one (E x [n\j = X 2 E x [n]), since J 3> SIE and SIE is in 
general positive. Moreover, due to the form of the func- 
tion f a , a symmetric fractional dissociation is always fa- 
vored. The most important result of Eq. (|62p is however 
the fact that the derealization error (or dissociation er- 
ror in ) has a twofold origin and only the second term 
on the right hand side depends on the self-interaction er- 
ror for the hydrogen atom. The remaining contribution 
traces back only to the wrong behavior of the functional 
under the homogeneous scaling. Note also that the first 
term might be the dominant one since the SIE is typically 
one order of magnitude smaller than the Coulomb energy 
in the H atom (see Table [I] for the value of SIE[nH\ for 
some popular functionals). 

This result shows the importance of the homogeneous 
density scaling in the development of accurate XC func- 
tionals. In fact, several functionals were constructed to 
minimize the self-interaction error of the hydrogen atom 
[94| , but no semilocal approximation exists that provides 
the correct scaling behavior under uniform density scal- 
ing. 

D. Fractional scaling (/? = 1) 

Let consider the disintegration of the hydrogen atom 
into two neutral atoms having a fractional number of elec- 
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trons and corresponding fractional nuclear charge. The 
electron density of the hydrogen atom with a fractional 
number of electrons q, must be considered as the ensem- 
ble density given by the superposition of the density of 
an hydrogen atom with one electron and nuclear charge 
q, with weight q, and that of one hydrogen atom with no 
electrons and nuclear charge q, with weight zero (this lat- 
ter thus does not contribute). Recalling that for the hy- 
drogen atom, a scaling of the nuclear charge corresponds 
to an uniform scaling of the density, we find that 



n,(r) = q 4 n 1 (qr) . 



(63) 



Generalizing the result of Eq. (|63l) , we define the frac- 
tional scaling as the scaling obtained from the family of 
relations given in Eq. (Q| when f3 = 1. Under this scaling 
the density and the various density parameters behave as 



n\(r) = A n(Ar) 
s A (r) = \~h (Ar) 



r sA (r) = A- 4 / 3 r s (Ar) 
q x {r) = \-iq{\r) , 
« A (r) = A 5/9 w(Ar) . 



(64) 
(65) 
(66) 



In the limit A — > oo the system is scaled towards the 
high-density limit (r s — > 0) and the universal functional 
is well approximated as 



F[n x ] 



-A 7 / 3 (E, 
-X 5/3 E. 



LDA r 



T. 



GB4r 



GE2\ 



(67) 



Note however that in this limit the reduced gradients for 
the correlation are not small (and thus the correlation en- 
ergy is not important). Eq. (|67|) shows a particular clas- 
sical behavior of the electrons: their kinetic energies be- 
come dominant over the classical Coulomb and exchange 
energies. Thus, in the limit A — > oo, the electronic system 
shows similarities with a gas of non-interacting particles, 
with non-uniform density ri\(x) — A 4 n(Ar). 

In the opposite limit (A < 1), the fractional scaling 
describes, as discussed above, the scaling towards a frac- 
tional atom. In fact, in this case the system is correctly 
scaled towards the low-density limit with the reduced 
gradient and Laplacian for the kinetic and exchange en- 
ergy becoming large. To strengthen the significance of 
such a scaling we consider its application to the disinte- 
gration of an atom, as introduced above. For simplicity 
we consider the disintegration of a hydrogen atom into 
two neutral atoms having a fractional number of electrons 
q and 1 — q, respectively, and corresponding fractional nu- 
clear charge. We define the XC disintegration energy as 



M(q) = E xc (l) - E xc (q) - E xc (l - q), 



(68) 



where E xc (q) denotes the XC energy of the H atom with 
fractional electron number and nuclear charge q. Note 
that correlation will only play a role when approximate 
non-self-interaction-free DFT functionals are considered 
(M DFT ), while only exchange will contribute in the com- 
putation of M exact . This quantity is very important since 
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FIG. 3: (Color online) Gaussian density XC error versus dis- 
integration error A (upper panel) and disintegration error A 
versus mean absolute errors of the AE6 test (lower panel), for 
several functionals (LDA [J], PBE [HJ, revPBE [§j|, APBE 
[lU, PBEint HU, and PBEsol @. The dashed lines are lin- 
ear interpolations of the data. 



the accuracy of any GGA functional in computing the 
values of M{q) is directly related to its ability to predict 
good atomization energies, because both processes (disin- 
tegration and atomization) preserve the total number of 
electrons. This fact is clearly shown in Fig. [3]where it can 
be noted the linear relation between the mean absolute 
error (MAE) on the computation of atomization energies 
of organic molecules (we considered here the AE6 test 
95], which is representative for organic molecule atom- 
ization energies) and the disintegration error A defined 
as 



A = 



M DFT (q) - M exact (q)] dq 



(69) 



for several representative XC functionals. Of course, a 
linear relation is obtained as well with the errors on the 
XC energy of a Gaussian one-electron density, which was 
shown to be a model system for atomization energies [35], 
£3. 

The disintegration error can be written alternatively 



E 



DFT, 



l)-EZ ac \l)-2 [E»"( q )-E: 



'(g)] dq. 
(70) 



For the exact energy we have E x xact (q) = E x xact (q) = 
q 3 E x xact (l), while for the approximate functionals 
we found numerically (see section |V| E xc FT (q) = 
q a E° FT (l) with a -2.3. Substituting into Eq. (JTOJ) we 
obtain 



A - 0.39K?/ T (1) 



_ E e T ct {1) 



(71) 
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Equation (JTTJ) shows that: (i) the hydrogen disintegra- 
tion problem is not related to the hydrogen one-electron 
density, as might appear at a first sight. Indeed, for a 
GGA which is exact for the hydrogen atom we would 
have A ~ 22 kcal/mol. In fact, Fig. [3] shows that the 
hydrogen atom disintegration is rather linearly related 
to the Gaussian one-electron density, (ii) The error A is 
determined basically from the behavior of the functional 
for q=l. (iii) To minimize the disintegration error (and 
thus yield good atomization energies), E E C FT must be 
about 10-20% larger than E^i act for the hydrogen atom; 
this condition is well met by GGA functionals accurate 
for atomization energies (APBE [5l|, revPBE J3Sf]), and 
is also well satisfied by the recently proposed zvPBEsol, 
and zvPBEint [35]. However, this result shows that in 
fact all these functionals base their performace on an er- 
ror cancellation effect which is made inevitable by their 
inability to respect the fractional scaling. 

Thus, the fractional scaling appears as an important 
exact constraint that XC functionals should respect. At 
present all popular approximations to the XC energy are 
only designed to respect the uniform scaling and are un- 
able to fulfill the fractional scaling relations. 

E. Strong-correlation scaling (/3 = —1) 

For /3 < —1/3, and in the semiclassical limit A — > 
oo, E LDA [n\] = X 0+4 / 3 E LDA [n] starts to dominate 
over the Thomas- Fermi kinetic energy T^ D [n\] = 
^2/3+5/3 rpLDA^ n ^ anc j ^ nus system may reach a 

strongly-correlated limit. Let us briefly consider the case 
p=-l. 

With this choice the scaling of the density and the 
reduced gradients are 

n A (r) - A- 2 n(A- x r) , r sA (r) = A 2 / 3 r s (A- 1 r),(72) 
s A (r) = X-^siX-h) , qx(r) = A" 2 / 3 ?^ -1 ') (73) 
i A (r) = A- 2 / 3 t(A- x r) , v x (r) = X^^v^r) (74) 

and thus, under this scaling with A — > oo, the system is 
in a slowly-varying, low-density limit. Thus, our approx- 
imation for the universal functional becomes 

F[n x ] « XJ[n]+X^{E^ DA [n]+E LDA [n]) 
A -i/3 (T s LDA [n] + E^ E2 [n}) + 
+\- 1 T° E2 [n] + X- 5 ^ E4 [n] , (75) 

which shows that the kinetic energy of the electrons, 
is much smaller than the Coulomb interaction, and the 
system resembles the features of a Wigner crystal, be- 
ing in a strongly-correlated limit. Note that in this 
limit, the LDA correlation energy scales as E LDA [n\] = 
X 1 / 3 E LDA [n], being as important as the exchange part. 
We recall that such a Wigner crystal is well des crib ed by 
a semilocal functional (see Eq. (23) of Ref. jlOl| ) de- 
rived from the point-charge-and-continuum (PC) model 



TABLE II: Effective scaling order s e ff (see Eq. ((76)) for sev- 
eral semilocal functionals at ft — and (8 = 1. The last line 
reports the reference value for the exact exchange functional. 
The second column reports the reference to the appropriate 
literature for each functional. 



Functional 


Ref. 


/9 = 


p = l 




Exchange-only functionals 




LDAx 


[86, 87] 


1.333 


2.666 


B88 


[91] 


1.254 


2.206 


PBEx 


f59] 


1.266 


2.323 


APBEx 


[5LJ 


1.262 


2.316 


revPBEx 


[88] 


1.252 


2.252 


PBEsolx 


[89] 


1.280 


2.356 


PBEintx 


[90] 


1.272 


2.329 


TPSSx 


[99] 


1.297 


2.399 


revTPSSx 


[100] 


1.309 


2.439 


Exact 




2.000 


3.000 




Exchange-correlation functionals 




t n* 
IjJJA 


[86, 87, 92] 


1.317 


2.01 ( 


PBE 


[59] 


1.272 


2.342 


APBE 


[51] 


1.268 


2.333 


revPBE 


[88] 


1.258 


2.271 


PBEsol 


[89] 


1.287 


2.378 


zPBEsol 


[89, 93] 


1.268 


2.316 


PBEint 


[90] 


1.279 


2.351 


zPBEint 


[90, 93] 


1.261 


2.289 


Exact 




2.000 


3.000 



[lOl| . Note that the PC model was inc orpor ated in high- 
level methods (e.g. ISI method of Ref. jl02| L and recent 
excellent work has been d one for fur ther development of 
the PC model, see Refs. flOl Hoi, as well for a DFT 
of strongly-correlated systems [l05l Il06| . Thus, further 
study of this scaling can be important. 

V. EFFECTIVE SCALING FOR SEMILOCAL 
FUNCTIONALS 

In this section we consider an assessment of semilo- 
cal XC density functionals for various scalings discussed 
above. In previous sections we discussed the scaling prop- 
erties of several exact energy functionals under the scal- 
ing transformations of the type of Eq. Q and we proved 
the utility of such scaling relations in various contexts. 
At the same time, we noted that these scaling relations 
are not respected by approximate functionals even at the 
LDA or second-order gradient-corrected level (the von 
Weizsacker kinetic energy functional being one excep- 
tion) . The situation is even worst for generalized gradient 
approximations which indeed do not have a well defined 
scaling behavior under the scaling transformations of Eq. 
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Nevertheless, it can be seen that most XC functionals 
in fact display an effective scaling E xc [n\] = X a E xc [n}. 
Solving this equation for the parameter a, we find a — 
[\n(\E xc [n x }\) - ln(|£ xc [n]|)]/m(A). Therefore, it is con- 
ceivable to define an effective scaling order for a generic 
functional E xc as 

f 1 \n(\E xc [n H }\) -ln(\E xc [n HX }\) ^ 
■seff = ^ W} dX . (76) 

Because our interest for the scaling of general exchange (- 
correlation) functionals in this work was motivated by 
the homogeneous and fractional scaling, we restricted our 
definition to use of hydrogen density nu — exp(— 2r)/ir 
and the interval A 6 (0 : 1), where we can compare with 
exact results. 

The effective scaling order provides a measure for the 
scaling behavior of different functionals, resembling in 
this respect the effective homogeneity of XC functionals 
[96Tj98j . Of course, for functionals having a well defined 
scaling behavior (e.g. LDAx), the effective scaling or- 
der will coincide with the analytic scaling exponent. For 
other functionals it will provide a measure of the effective 
scaling behavior, so that the deviations of the effective 
scaling order from the true value could give an estima- 
tion of the accuracy of the functional to fulfill the scaling 
relation. We note that for the latter cases the integrand 
of Eq. (|76p was always found to be almost constant over 
the entire integration interval (except very close to the 
boundaries; note however that the function is integrable 
over the given range) , showing the robustness of our def- 
inition. 

In Table [TTI we report the values of s e // for several ex- 
change and exchange-correlation functionals at /3 = and 
/3 = 1 . An inspection of the data shows that all the func- 
tionals perform similarly and quite differently from the 
exact reference, that in this case is the exact Kohn-Sham 
exchange, since for n — uh and A < 1 there is no correla- 
tion. Remarkably, the best scaling behavior is found for 
LDA exchange, while slightly worst results are obtained 
for GGA functionals. The use o f meta-GGA functionals, 
as TPSS HI or revTPSS [ijxj, which are constructed 
taking into account the physics of one-electron systems, 
is found finally to bring a slight improvement in the ef- 
fective scaling behavior. Moreover, at the GGA level 
the addition of approximate correlation seems to bring 
some small improvement in the scaling behavior, in line 
with the fact that indeed semilocal DFT functionals are 
not really exchange or correlation functionals but rather 
rely on an heavy error compensation between the two. 
We recall instead that meta-GGA functionals are one- 
electron-self-correlation free. 



VI. CONCLUSIONS 

In summary, we have investigated the scaling with vari- 
able particle number of the form of Eq. n\{r) = 



\ 3,3+1 n (\Pr). For such scaling transformations we pro- 
vided a formal definition within the ensemble formalism 
of DFT and studied the basic features, also in relation 
to the scaling properties of different important density 
functionals. 

The density scalings defined in Eq. @, spans an im- 
pressive set of physical properties: in the limit of large A 
they are crucial for semiclassical theory of many-electron 
systems (e.g. Thomas-Fermi scalings is related to atoms, 
Uniform-electron-gas scaling is related to metallic clus- 
ters, strong-correlation scaling is related to Wigner crys- 
tals), whereas in the limit of small A they are connected 
with the physics of small systems with fractional particle 
number, and to self-interaction errors. 

The here proposed uniform-electron-gas scaling ((3 — 

— 1/3) is the right basic concept for jellium clusters. Sim- 
ple scaling manipulations showed that the curvature cor- 
rections are in fact related to the second-order gradient 
expansion. On the other hand, the surface corrections 
(described by the Airy gas model (76[), are quantum os- 
cillations terms. By analo gy with the recent work on 
the semiclassical atom [26l l49l l50j . a modified second- 
order gradient expansion (MGE2) can be constructed for 
jellium clusters, in order to recover the exact surface cor- 
rections. However, we expect that such a MGE2 will be 
very close to the regular GE2 that is accurate for surfaces 
of simple metals when the Kohn-Sham densities are used 
89]. * 

Moreover, the idea of MGE2, that can account for the 
principal quantum corrections, can be generalized for any 
P (in the limit of large A). Such a /3MGE2 will be very 
useful especially for the strongly-correlated scaling (/3 = 

— 1) where the LDA term is not exact (in the limit A — > 
oo ) d ue to the self-interaction problem in Wigner crystals 
[lOlj - For example, /3MGE2 for the exchange energy may 
have the form 

£f MGE2 M = jdvn e L x DA f({3)(l + M (/3) S 2 ), (77) 

where /(/3) and fJ,(/3) should be derived in further investi- 
gations. (For the Thomas-Fermi scaling f(/3 = 1/3) = 1 
and [j,(l3 = 1/3) = 0.26 [H; whereas the regular GE2 has 
/ = 1 and n = 0.12346). 

Finally we have shown the usefulness of the here pro- 
posed fractional scaling (/3 = 1) for the atomization en- 
ergies of molecules. Recently, it has been derived an 
atomization energy constraint (i.e. minimization of an 
entropy-like function for an ensemble one-electron den- 
sity models) (35l . |93| . This constraint was derived from 
an empirical observation relating errors in the model one- 
electron densities to errors in the atomization energies, 
of popular GGAs (see Fig. 2 of Ref. [93j); and from the 
physical explanation that one-electron densities are sim- 
ple models for simple bonding regions, where iso-orbital 
regime can be significant, (see Fig. 1 of Ref. [93j, and 
the corresponding discussion). Using the fractional scal- 
ing (/3 = 1), we have better explained the significance of 
one-electron Gaussian model for atomization energy of 
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molecules and disintegration of the hydrogen atom (see 
Eq. (fTTj) and its related discussion). 

We recall that the semilocal exchange hole models sat- 
isfy the s um rule for systems with integer number of 
electrons |108| . but violate the exchange hole sum rule 
in case of fractional number of electrons, and thu s pre- 
dicting too- negative energies for such systems [l09j | . The 
here proposed effective scaling (see Section V) , measures 
in fact the functional accuracy for systems with fractional 
particle number, and can be used to develop (and test) 
new better approximations. 

Our work provides a deeper insight into the relevance 
of the scaling relations having the form defined in Eq. 

and highlights the importance of these scaling rela- 
tions in DFT. In particular, Eqs. (|3T j) - ([34| provide use- 
ful scaling relations with varying particle number that 
are important constraints in the construction of approx- 
imate exchange-correlation (or noninteracting kinetic) 
functionals. Moreover, for the A — > limit, Eqs. ([57)1 
and are exact constraints for the Kohn-Sham ex- 
change and kinetic energy functionals, respectively and 
the effective scaling order (Eq. ([75)) ) is a more general 
requirement, which is relevant for the SIE problem. Un- 
fortunately, no such explicit expressions exist instead for 
the A — >■ oo limit, where however the recovery of the 
scmiclassical atom physics was shown to be an impor- 
tant condition. 
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Appendix A: Scaling relations for Kohn-Sham 
kinetic and exchange energies for one electron 
systems 

In this appendix we consider briefly the special case of 
a fractional number of particles (with N < 1) and A < 1. 
In this case, using the formalism of Section|TT]we can de- 
rive, in line with Ref. 0, the general scaling properties for 
the non-interacting kinetic energy and the Kohn-Sham 
exchange functionals. Under coordinate scaling the den- 
sity operator is transformed into another valid density 
operator and because of the simple scaling behavior of 
the statistical weight we have f A (r) = Af (A^r). At this 



point we can write Tr 



A 2/3+l Tr 



TT 



where we 



used the fact that T(r/A' 3 ) = \ 2l3 f(r), and we note that, 
for any A, if Tr TT 
Hence, 



is a minimum so must be Tr 



T a [n x ] = X 2fi+1 T s [n} . 
In a similar way it can be proved that 
E x [n x ]=\?+ 2 E x [n] . 
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